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Abstract 

We present exact expressions for the eigenvalues and eigenvectors of 
the d dimensional Laplace operator in a cut Fock basis. 

Introduction 

The ultimate goal of studies involving Supersymmetric Yang-Mills Quantum 
Mechanics (SYMQM) is the solution of the system defined in D = 9 + 1 space- 
time dimensions with SU(N) gauge group for any N. The latter appears in 
such problems as the quantization of a supermembrane [1] [2] or the description 
of dynamics of DO branes, one of the kinds of objects present in M-theory [3] [4]. 
It should not be surprising that finding such solution is a difficult task. Therefore 
one tries to first investigate the lower dimensional versions of SYMQM. Among 
them, the D = 1 + 1 SYMQM is the simplest one [5] . It is described by position 
<f> m and momentum n m bosonic operators and, A m and A" 1 fcrmionic operators. 
The index m labels the color degrees of freedom, m — 1,...,T>, where V is 
the dimension of the adjoint representation of the gauge group. Although the 
D = 2 SYMQM Hamiltonian is free, its solutions are not trivial due to the 
singlet constraint which represents the dimensionally reduced Gauss law. In 
other words, the physical Hilbert space of D = 2 SYMQM is composed of states 
invariant under the gauge transformations. 

A convenient method for analyzing such systems was proposed few years 
ago by Wosiek [6] (for a recent introductory article see [7]). It was designed 
mainly as a fully non-perturbative, numerical method for deriving the spectra of 
SYMQM. One of its advantages is that it treats fermions and bosons on an equal 
footing which made the calculations in all fermionic sectors possible. Moreover, 
it can handle any gauge group in any dimension. The main idea is to calculate 
the matrix elements of the Hamiltonian operator in a Fock basis composed of 



*e-mail address: korcyl@th.if.uj.edu.pl 



states containing less than Nb quanta. Nb is called the cutoff. Such finite 
matrix is then diagonalized numerically in order to obtain the approximated 
eigenenergies. A suitable choice of the cutoff, preserving the gauge symmetry 
as well as the rotational symmetry, enables one to obtain many numerical results 
M BE] H5 • The cut Fock space approach turned out also to be a suitable tool 
for analytic investigations. 

In this short paper we analyze the D = 2 SYMQM models with the SO(d) 
gauge group. In the language of SYMQM they can be interpreted as an ex- 
tremely simplified c£-dimensional model with only one color or, equivalently, as 
a model with d colors transforming in the adjoint representation of the SO(d) 
group. Due to the simplicity of the SO(d) symmetry group we are able to de- 
rive analytic expressions for the cut spectrum and eigenstates. Consequently, 
we check that the infinite cutoff limit of the cut solutions is given by the well- 
know Bessel functions [H][I3], which are solutions of the free d-dimensional 
Schrodinger equation carrying zero angular momentum. We expect that such 
analytic solutions can be used as a starting point for the study of systems with 
several colors transforming with the SU(N) group. Moreover, once known for 
any d or N they could be instructive in the further study of physically interesting 
'large-d' or 'large-N' limits. 

We also discuss the physical interpretation of the cutoff imposed on the 
Fock basis. We show that it corresponds to a specific discretization of position 
and momentum spaces. Afterwards, we briefly analyze the scaling law, i.e. 
the dependence of the eigenenergies on the labeling index. Such law was first 
investigated by Wosiek and Trzetrzelewski [14] [15] in one spatial dimension and 
was used as a handle to correctly differentiate localized states from the non- 
localized ones within our numerical approach. In this note we discuss its validity 
in higher dimensional spaces. 

We summarize our results in the last section. 

The system 

Traditionally, in the D = 2 SYMQM models cp m and tt™, as well as, A m and 
A™ transform in the adjoint representation of the SU(N) gauge symmetry. In 
this paper we concentrate on the SO(d) gauge symmetry, and therefore, we 
assume <p m and A'™ to be in the adjoint representation of SO(d). The SYMQM 
Hamiltonian is supplemented by the singlet constraint, hence we search for 
solutions carrying zero angular momentum. The Hamiltonian reads 
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fulfilling the commutation relation, 

[a",a tm ] =S mn , (3) 

and rewrite H as 

H = a) m a m + -- -a m a m - -a^ m a) m , (4) 
2 2 2 w 

where the summation from 1 to d over repeated indices is understood from now 
on. 

The SO(d) group possess only one invariant symmetric tensor, namely the 
Kronecker delta. Hence, the algebra of invariant under rotations creation oper- 
ators reads 

[a m a m ,a) n a) n ] = 2d + ia^ n a n , 

[a tm a"\ a + V"] = 2o tn a tn , (5) 
[a p ,a tm a tm ] = 2a tp . 

The basis vectors in the singlet sector can be labeled by a single index n, 
and take the form 

|n) = i(at-at-)"|0>. (6) 

Notice that n is related to the half of the total number of quanta of each type 
contained in \n), Nb — 2n. For a given cutoff Nb only states \n < [^-\) arc 
present. There are [^"J + 1 such states, a new one appearing for Nb even. 

We calculate the normalization factor in eq. © by commuting consecutively 
all operators a q a q through (a' p a' p )' 1 , using eqs.©, so that they annihilate the 
Fock vacuum. This gives a recursive relation for A/"", 

(Kf - 2n(2n + d 2)«_ 1 ) 2 I «) 2 = 1, (7) 
which is solved by 

(<) 2 -4«n!£|±^. (8) 

Although we are mainly interested in the scalar sector, we can obtain results 
for higher representations without additional effort. As an example, we present 
some calculations in the vector sector. The basis states in this sector can be 
obtained out of the states given by eq.© by the action of an additional single 
creation operator, 



|n,i) = iatV m a tm ) n |0), i = l,...,d, (9) 



where 

(Af}) 2 = 4 n n!- 

r(f + 1 



l^) 2 = 4"n! r( j+" + 1) . (10) 
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States given by eq.© are orthogonal in n and i: (m,j\n, i) = S mn Sji(M^) . 

It is now straightforward to evaluate all elements of the Hamiltonian matrix 
in the singlet sector, 

(n\H\n) = 2n+ |, (n + l\H\n) = (n|# |n + 1)=-J(n + 1)(| + n), ( u ) 
and in the vector sector, 

(n,j|fl"|n,i> =%(- + 2ri+l), 

(n,i|F|n + l,i> = (n + l,j|2T|n,i) 

This done, we can calculate numerically the spectra. In order to obtain 
the physical results corresponding to the infinite cutoff limit we investigate the 
behavior of the eigenenergies with increasing cutoff for some finite values of Ng . 
Figure[T]contains such results for the 12 lowest eigenenergies in the singlet sector. 
Each curve falls slowly to zero as expected from the nonlocalized character of the 
eigenstates of the free Hamiltonian [14] [15] [16] . The dependence on the cutoff 
will be also discussed in the last section of this note. 
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Figure 1: Spectrum of the Laplace operator in d = 3 dimensions in the singlet 
sector. 
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Character method 

There exists an independent way to calculate the number of basis states trans- 
forming according to a given representation of the SO(d) group, which exploits 
the orthogonality of characters of irreducible representation^. We start with a 
formula for the number of singlets appearing in a product of Nb adjoint repre- 
sentations [T7] 

4b = / d PffoM (oi) 1 ■ X% B {ai), (13) 

where d^ so{d) (on) is the invariant Haar measure for the SO(d) group parame- 
terized by the set of angles {o;i} i _ 1 \d/2\ runnm § from to 2ir and X% B { a i) 
represents the symmetrized product of Nb characters of the adjoint representa- 
tion of the SO(d) group and is equal to [T5] 

N B 

xk(«0= E IlrT^(^W)"' ( 14 ) 

Y, k ki k =N B k=l 

1 in eq. (fT3")) stands for the trivial character of the invariant representation of 
SO(d). It can be replaced by the character of any other representation R, in 
which case, d% will be equal to the number of times R appears in the product 
of Nb adjoint representations. For d even, d — 2M with M integer, the Haar 
measure for SO(d) is given by [19] 



2 {M-xy , , 2 

d nso W = f M M! 11 [cos(a k ) -cos(oy)J , (15) 



2 (M-1) 

r 

l<j<k<M 

and for d odd, d = 2M + 1, 



2_M 2 2 M 

d nso W = ^M^j II (cos(a fe ) -cos(aj)) Y[sm 2 {^-). (16) 

l<j<k<M m=l 

The characters of the adjoint representation of SO(d) are given by, for d = 2M 

M 

X d K) = 2 ^ cos(a m ), (17) 



and for d = 2M + 1, 



a/ 



^(a,) = 1 + 2 ^ cos(a m ). (18) 



m — 1 



Explicit evaluation of the integrals confirms that eqs. ([6]) and (|9]) form indeed 
a complete basis of the Hilbert space in the corresponding sectors. 



The method was suggested by R. Janik 
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Analytic results 



In this section we derive the set of eigenvalues and eigenvectors of the Hamilto- 
nian operator in the singlet and vector sectors for any finite cutoff Nb- 

At cutoff equal to Nb, the basis size in the singlet sector is n max — |_ 2 J + 1- 
Let us denote by In max (A) the characteristic polynomial of the Hamiltonian 
matrix at this cutoff. Then, I\ = (0\H — A|0) = | — A. Due to the tridiagonal 
nature of the Hamiltonian matrix we can write a recursive relation for I n (A), 

4(A) = (2n - 2 + | - A) J n _i(A) - ((n - l)(n - 2 + ^))/„_ a (A). (19) 

By changing the variables as I n (X) — ► ra!J„(A) we get a recursive equation which 
is solved by generalized Laguerre polynomials, i.e. 

I„(A) =i| _1 (A). (20) 

These well-known polynomials are defined as the solutions of the differential 
equation 

xy" + (a + 1 — x)y + ny = 0, 
and the orthogonality relation 



L^ n (x)L^(x)x a e x dx = S mn . 

Therefore the eigenvalues of the Hamiltonian operator in the singlet sector in 
d dimensions at cutoff Nb are given by the zeros of the Laguerre polynomial 

L\ B (A). By taking d = 1 we recover the results of [TS]. 

An alternative method to derive this result is to consider a general state \E) 
from the singlet sector and expand it as 

\E) = £ a n (£)|n>. (21) 

The eigenequation of H translates into a recursive relation for the a n {E) coef- 
ficients. After extracting an irrelevant constant factor, a n {E) — > 2~ n a n (E), the 
recursion reads 

- (2n+ - - E)a n {E) + (n + l)(n + ^)a„+i(S) = 0. (22) 
Eq.p2) is solved by 

a n {E) = a T{t)^^ y n = l,..., L ^j ) (23) 

--i 

together with the quantization condition L\ E ^(E) = 0, which is equivalent 
to eq.O- 
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Similarly, in the vector sector the recursive relation for the determinant of 
the Hamiltonian matrix can be solved by 

InW = iJ+i(A). (24) 

Thus, the spectrum in the vector sector for the cutoff Nb is given by the zeros 
of the generalized Laguerre polynomials with the index shifted by 1 comparing 
to the singlet case, 

& TVr 
R\E\) = E n \E l n ), with E n such that L* =0, n = 1, . . . , + 1 

where \E l n ) are the eigenstates of H with the vector index i. As a straightforward 
consequence of the vector nature of the states each eigenvalue is d times 
degenerate. The corresponding eigenvectors are given as an expansion in the 
basis 

\E l ) = £ a nti (E)\n,i). (25) 

n=0 

with 

r( d \ L ^ E ) i r^l loa\ 
a n AE) = anV(-) -r 1 , n=l,...,[ 1. (26) 

,V 2 ; r(n+| + l) ' 2 1 v ' 

Hence, we obtained simple formulae giving the spectrum of the cut Hamilto- 
nian operator as well as its eigenstates. The dependence on the cutoff is explicit 
and it is straightforward to perform the infinite cutoff limit. 

One can easily generalize the above expressions for the eigenvalues and eigen- 
vectors transforming under any irreducible representation of the SO(d) group. 

Reconstruction of radial waves 

In this section we confirm the results obtained above by comparing their infi- 
nite cutoff limit with the well-known continuous results. Before we present the 
detailed calculations, we make two observations. 

The rescaled Hamiltonian of the quantum harmonic oscillator is given by 

H osc =p 2 + x 2 . (27) 

It is evident that H osc is invariant under the transformation x o p. Therefore, 
the functional form of the wave functions in the momentum and position repre- 
sentations is the same up to a constant multiplicative factor. In fact, both are 
related by the d dimensional Fourier transform. 

Next, the amplitude a n {E) can be interpreted as the n-th wave function 
of the (i-dimcnsional harmonic oscillator in the momentum representation. To 
see this, let as denote the plane wave with a wave-vector k, k = V~E, which is 
the eigenstate of the free Hamiltonian H, by |/t), and the n-th state of the d- 
dimensional harmonic oscillator carrying zero angular momentum by \n). Then, 
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a n (E) is simply the scalar product of the free wave \k) and the harmonic oscil- 
lator eigenfunction |n) , 

a n {E) =a n (K 2 ) = ( K \n). (28) 

Thus, 

ip n ( K ) = (K\n) = /(K 2 )i! _1 (K 2 ), (29) 

cj )n (r) = (r\n)=g(r 2 )Lt 1 (r 2 ), (30) 

where f(n 2 ) and g(r 2 ) are functions which do not depend on n and could not be 
determined from the recursion relations. We can check the above conclusions 
by explicitly solving the radial d dimensional harmonic oscillator eigenequation 
in the position representation, 

ld 2 (j) n (r) d-ld(j) n (r) d 1 2 

+ ~r dV~ + {2n + 2 )0 " W ~ 2 r Mr) = °' (31) 

The solutions to this equation are 

Mr) = Cl e-4r 2 -\ Fl (l - \ - |, 2 - * r 2 ) + c^^- |, * r 2 ). (32) 

The first term of the solution is singular at r = for d > 2, so we set c\ = 0. 
By rewriting <p n (r) in a more familiar form we get 

Mr) =c 2 e-^4~V), (33) 

and similarly 

= c 2 e-T J L|" 1 (K 2 ), (34) 
where the constant c 2 can be fixed by the normalization condition, 



f'2 



(35) 



\ 7rfr(n+f)' 



At this point we can calculate the wave functions in the position representation 
of the infinite cutoff limit of solutions eq. (|23|) . We write 



(r|«) =^(r|n)(n|«) = ^ ^nW^W 



n=0 



K 2 +r 2 «-—« f— T)™rjl a i * 1 

_ 1 + 2 J 



with all unimportant constant factors gathered in Af and the (—1)™ factor come- 
ing from the definition of Fourier transform. Using the known relation for asso- 
ciated Laguerre polynomials [13] we eventually get 

(r\K) = l -N- L—J^-Kr). (37) 

z (—Kr)2 



8 



This result indeed coincides with the solutions of the radial part of the Laplace 
equation, 

- A*(r) - Eftr) = *M + *=1 « + ^(r) = 0. (38) 
ar z r ar 

given by the well-known Bessel functions, 



Cir 1 



i Ji_ x {-\/Er) +c 2 r 1 -iY< L _ l {-\ r Er). (39) 



Requiring an analytic behavior of the solution at r = fixes C2 = 0, and 
after setting c\ = ^A/Vc 1- ^ and E — n 2 , we recover the solutions of eq. ([37]) . 
Therefore, our solutions obtained for a finite Nb become, in the infinite cutoff 
limit, the known Bessel functions. 



Physical interpretation of the cutoff 

Physically, the introduction of a cutoff limiting the maximal number of quanta 
contained in the basis states is equivalent to a specific discretization of the posi- 
tion as well as momentum variables. This can be demonstrated by considering 
the eigenvalues of the cut matrix representations of those operators. They cor- 
respond to possible outcomes of a measurement process. Obviously, in the cut 
basis the matrices of the x and p operators are finite and have as many eigen- 
values as there are states in the cut basis. Hence, position as well as momentum 
spaces are discrete. 

In order to be more precise, let us consider gauge invariant position and 
momentum operators of the form X 2 = x m x m and P 2 = p m p m . We have 
shown in the preceding sections that their spectra consist of the zeros of Laguerre 
polynomials. We denote their eigenvalues by X 2 and P 2 , respectively, 

X?: Lr\xf)=0, P 2 : Lr\p?) = 0, i = l,...,n. 

A well-known upper bound for the largest zero of L"(x) is, provided \a\ > j, 

a > -I, 

x n < (V4ra + 2a + 2 - 7(4n + 2a + 2)~^) , 

where 7 = 6~iii and i\ is the smallest zero of the Airy function [20]. Hence, 
for a cutoff Nb the eigenvalues Xi and Pi lie in the interval I, 

X=(-^l + 3, f^l+3). 

The interval I of length ^/Nb contains roughly \Nb such eigenvalues. There- 
fore, one concludes that and Pi tend to cover the real axis in the limit of 
infinite cutoff. We call this limit also a continuum limit. 
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Continuum limit 



Up to now, our results suggest that the continuum limit of the n-th eigenenergy 
vanishes, 

lim E n {N B ) = 0, (40) 

N B — >oo 

However this is in contradiction with the general expectation that the spectrum 
of the kinetic energy operator is continuous, i.e. 

lim E n (N B ) = \p 2 , pel, (41) 

In order to avoid the situation of ea. (|40p one must assume that the limit in 
eq. (|4*Tj) is taken with, 

n = n(N B ), (42) 

which we call scaling law, referring to a similar procedure in lattice field theory. 
This phenomenon was investigated in the case of one-dimensional momentum 
operators by Trzetrzelewski [15]. He was able to derive the scaling properties 
of the eigenvalues of the momentum operator as well as of the kinetic energy 
operator. We report the latter in table [TJ The four cases shown correspond 
to different parities of the corresponding eigenstates. Such complication is not 
present for d > 1 since the basis states are always even, i.e. build of powers of 
at 2 . 





N B odd 


N B even 


n odd 


E„{N B ) w 2 2Nb+3 




n even 


E n (N B ) 2 2Nb+5 


E n (N B ) ~ 2 2N B +3 



Table 1: Scaling laws for the eigenenergies for even and odd N B and n. 

It is now straightforward to check the validity of the scaling law in d > 1 
dimensions. The eigenenergies of the kinetic energy operator being the zeros 
of appropriate Laguerre polynomials can be approximated at large N B by the 
formula [12] [13] 

H4^ 1 W)„- i J^ + ..-, (43) 

Hence, we see that the simple scaling E n ~ n 2 is correct as long as the following 
approximation is true 

i$_i,„»7i(<0n + '»(<9- ( 44 ) 
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Figure 2: Scaling laws for d = 1 and d = 150. 

Obviously, ([44)) will not hold for high d. Indeed, Figure [2] shows the depen- 
dence of the eigenenergies E n (Ns) on the labeling number n for a given cutoff 
Nb = 200. Distinct curves correspond to different d, namely d = 1 for the lower 
one and d = 150 for the upper one. Definitely, in the latter case the dependence 
is no longer linear. Some particular values of the coefficients 71 and 72 are given 
in table [2] 



d 


7i (d) 


Std. error 


12(d) 


Std. error 


1 


3.14159 


< io- llJ 


-1.5708 


< 10" 1U 


3 


3.14159 


< 10- 10 


0.0 


< 10- 10 


9 


3.142 


4.4 x 10~ 5 


4.6192 


0.0077 



Table 2: Fitted values of parameters 71(d) and 72(d) for d = 1,3, 9. 



Conclusions 

In this note we used the cut Fock space approach to investigate the system 
which can be interpreted as a free quantum particle in d dimensions. Although 
solutions to such problem are well known in the continuum (Nb = 00), they 
were unknown for finite Nb- 

Our approach was guided by the numerical algorithm proposed for solving 
supersymmetric Yang-Mills quantum mechanics in the Hamiltonian formulation. 
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The study of this particular system was motivated by a similar problem with 
the SO(d) gauge group replaced by the SU(N) group [21] [22]. 

We started by describing the construction of basis states with the use of 
gauge invariant creation operators. Then, we explicitly derived the Hamiltonian 
matrix. Due to its simple form it was possible to obtain exact expressions for 
the eigenvalues and eigenvectors. The formulae are parameterized by the cutoff 
and hence, are valid for any finite cutoff. Namely, the spectrum at finite cutoff 
is given by the zeros of an appropriate Laguerre polynomial, whereas in the 
continuum it corresponds to the positive real axis. 

In order to confirm the correctness of our solutions, we calculated their wave 
functions in the position representation in the continuum limit. Indeed, it was 
possible to show that such wave functions are equivalent to the well known 
solutions given in terms of Bessel functions. 

Next, we discussed the physical interpretation of the cutoff. We argued that 
a finite cutoff corresponds to a discretization of both position and momentum 
space. In the infinite cutoff limit continuum results are recovered. 

Finally, we observed that the scaling law proved for the one dimensional 
kinetic energy operator holds only approximately in spaces with higher number 
of spatial dimensions. 
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